Article 


Evidence for the earliest structural use of 
wood at least 476,000 years ago 


https://doi.org/10.1038/s41586-023-06557-9 


Received: 21 January 2023 


Accepted: 21 August 2023 


Published online: 20 September 2023 


Open access 


® Check for updates 


L. Barham'™, G. A. T. Duller?, I. Candy%, C. Scott’, C. R. Cartwright®, J. R. Peterson’, 
C. Kabukcu**, M. S. Chapot?, F. Melia‘, V. Rots’, N. George’, N. Taipale’, P. Gethin’ & 
P. Nkombweê 


Wood artefacts rarely survive from the Early Stone Age since they require exceptional 
conditions for preservation; consequently, we have limited information about 

when and how hominins used this basic raw material’. We report here on the earliest 
evidence for structural use of wood in the archaeological record. Waterlogged 
deposits at the archaeological site of Kalambo Falls, Zambia, dated by luminescence 
to at least 476 + 23 kyr ago (ka), preserved two interlocking logs joined transversely 
by an intentionally cut notch. This construction has no known parallels in the African 
or Eurasian Palaeolithic. The earliest known wood artefact is a fragment of polished 
plank from the Acheulean site of Gesher Benot Ya’aqov, Israel, more than 780 ka 
(refs. 2,3). Wooden tools for foraging and hunting appear 400 ka in Europe* ®, China? 
and possibly Africa’®. At Kalambo we also recovered four wood tools from 390 kato 
324 ka, including a wedge, digging stick, cut log and notched branch. The finds show 
an unexpected early diversity of forms and the capacity to shape tree trunks into large 
combined structures. These new data not only extend the age range of woodworking 
in Africa but expand our understanding of the technical cognition of early hominins”, 


forcing re-examination of the use of trees in the history of technology”. 


Inthe African context, indirect evidence for woodworking comes from 
use-wear traces and residues on Early Pleistocene stone tools in East 
Africa (Oldowan, Acheulean)“*™. Actual wood objects are found in 
Mid-Pleistocene waterlogged deposits in southern Africa with Acheu- 
lean and Middle Stone Age tools. At Kalambo Falls, wood was recovered 
from Acheulean horizons inthe 1950s-1960s (sites A and B), but tapho- 
nomic processes removed evidence of intentional shaping from most 
pieces”. A wood chip and three objects with transverse notches raised 
the possibility of intentional modification. Attempts to date the wood 
gave minimum ages”, At Amanzi Springs, South Africa, a single stick 
with a possible chop mark was reported from waterlogged Acheulean 
deposits excavated in the 1960s”. The deposits were radiometrically 
dated (approximately 404-390 kyr), with wood found in recent excava- 
tions but without evidence of modification”. The earliest clearly modi- 
fied wood object, collected in 1952 from spring deposits at Florisbad, 
South Africa, was associated with Middle Stone Age tools and hominin 
remains (Homo helmei)”. The object’s tip shows cutmarks and fine 
striations”, but its location relative to dated deposits is uncertain”. 
Excavations at Kalambo Falls in 2019 recovered five modified wood 
objects at site BLB (Fig. 1b) from four areas (BLB2, BLB3, BLB4 and BLBS) 
in sediments above and below river level (Fig. 1c). A sixth object, from 
BLB3, showed no evidence of modification. Two objects were associated 
with Acheulean artefacts below the river (BLB3 and BLBS); three from 
contexts above river level, without stone tools (BLB2 and BLB4) (Fig. 2). 


The Quaternary sequence is a 9-m-deep exposure above the Kalambo 
River (BLB1is a geological section). Sediments are fluvial sands and 
gravels with occasional, discontinuous beds of fine sands, silts and 
clays with wood preserved in the lowermost 2 m (Methods). A perma- 
nently elevated water table has preserved wood and plant remains 
(Supplementary Information Section 1). The depositional sequence 
is typical of a high- to moderate-energy sandbed river that underwent 
lateral migration. The sands are dominated by a lower unit of horizontal 
bedding and an upper unit of planar/trough cross-bedding. Upper and 
lower sand units are separated by fine sands, silts and clays with plant 
material deposited in still water after the river migrated/avulsed else- 
where in the floodplain. Wood is deposited in this environment either 
through anthropogenic emplacement, or naturally transported inthe 
flow, and snagged on sand bedforms (Supplementary Information 
Section 1). 

Dating is based on16 sand samples collected for luminescence analy- 
ses from deposits bracketing key finds including those containing 
wood. Younger samples are dated using single-grain quartz optically 
stimulated luminescence (OSL) and older samples by postinfrared 
infrared stimulated luminescence (pIR IRSL) from potassium-rich 
feldspars (Methods and Supplementary Information Section 2). The 
pIRIRSL approach used extensively in recent years”*”° does not suffer 
the problems that can generate large uncertainties associated with 
thermally transferred OSL (TT-OSL), as seen at Site C North (Fig. 1b)”°. 
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Fig. 1| Location of Kalambo Falls archaeological site and excavated areas. 

a, Site location in south-central Africa. b, Course of the Kalambo River (in outline) 
from around 1956 to 2006 in relation to previous excavations at sites A,B, C, D 
and C North. Site BLB (2019) and excavation units BLB1, BLB2, BLB3, BLB4 and 
BLBS are located along the current main channel (blue). c, Cross-section of the 
2019 excavation units showing the location of 16 luminescence dating samples 
(KFO1-KF17, dark blue circles, uncertainties (+) shown at 1- o) by unit. Unit BLB1 
isa geological section of the full cliff exposure from ground surface to below 
water level. The three colour bands indicate clusters of pIRIRSL ages grouped 


The new ages (Extended Data Table 1) are in stratigraphic order (10) 
(the only exception being KF10 and KF11), forming three clusters based 
on mean ages (and error of the mean) (Fig. 1c). The earliest cluster 
(476 + 23 kyr) incorporates wood in deposits below river level (BLB3, 
BLBS). The intermediate cluster (390 + 25 kyr) encompasses one wood 
tool in BLB2 found above river level. The upper cluster (324 + 15 kyr) 
brackets two wood objects in BLB4, above river level. No wood was 
found higher up the sequence. 

The wood was sampled for identification to species level (Methods 
and Supplementary Information Section 3)” and for radiocarbon dat- 
ing, to assess for an intrusive origin; all dates were infinite (more than 
50 kyr ago (ka), T. Higham personal communication). Infrared spec- 
troscopy shows partial mineralization of the wood (silica) (Methods 
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by mean ages and standard error of the mean. The earliest wood objects 
(BLB5, BLB3) are in the lower green band witha mean age of 476 + 23 kyr. The 
blue band has a mean age of 390 + 25 kyr and incorporates the wood object in 
BLB2. The overlying yellow band has a mean age of 324 + 15 kyr and incorporates 
wood objects in BLB4. Red diamonds indicate modified wood objects. Inset map 
data ina: Google Maps, Google 2021, INEGI; the map of Kalambo Falls was drawn 
using Arc-GIS Open Street Map. Map data in b: Google Imagery 2021 CNES/ 
Airbus, Maxar Technologies; the 1960s river course, with site locations, was 
redrawn from figure 7.2 of ref. 17 with permission from Cambridge Univ. Press. 


and Supplementary Information Section 5). Surface modifications 
were identified on photographically generated models (Methods). The 
need to keep the wood wet made standard photography, structured 
light scanning and microscopy impractical because of high surface 
reflectance and pooling of water in key features, obstacles overcome 
by photogrammetry and reflectance transformation imaging (RTI) on 
the submerged wood. Additional photographs were taken after brief 
episodes of drying. 

Manufacture marks and use history are described using standardized 
terminology’ with interpretations supported by replication experi- 
ments (Methods and Supplementary Information Sections 6 and 7). 
Functional interpretations are drawn from wood artefacts in water- 
logged Holocene sites in Zambia”, the UK”? and ethnographic sources 


Fig. 2 | Modified wood tools from site BLB, Kalambo Falls, 2019. a, BLB5 
structural element (object 1033). b, BLB3 ‘wedge’ (object 660). c, BLB2 ‘digging 
stick’ (object 219). d, BLB4 cut log. e, BLB4, tapered piece with single chop-mark. 
Scale bars, 10cm. 


(Methods and Supplementary Information Section 8)”. Individual 
pieces are described by excavation block from oldest to youngest. 

In BLBS, one modified object was recovered from basal sands contain- 
ing flake tools, cleavers and handaxes. Object 1033 is alog (Combretum 
zeyheri) 141.3 cm long x 25.6 cm, with tapering ends in three parts, 
which overlies a larger treetrunk at a 75° angle (Fig. 3a,b). The area of 
overlap is a U-shaped notch 13.2 cm long by 11.4 cm wide, transverse 
to the long axis (Fig. 4). The underlying trunk, also modified, and left 
in situ, passes through the notch (Fig. 3a). 

Bark (on parts 1 and 2) and a knot (part 1) occur near sapwood on 
object 1033. The notch (part 1) exposes tangential and radial longi- 
tudinal sections, with sapwood missing. Evidence of chopping and 
scraping occurs on the upper surface of the notch (Fig. 4, areas a-d, 
and Extended Data Fig. 1) and on the lower surface (Fig. 4, areas e-k, 
and Extended Data Fig. 1). Areaa has two transverse intersecting sets of 
scraping marks with superimposed parallel linear striations (maximum 
length 12.5 mm with V-shaped cross-sections, maximum width 3.2 mm). 
Area b has two linear transverse entry facets (18.8 mm long x 1.6 mm 
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wide). Area c is a single facet with clear entry and stop marks (24 mm 
long x 1.4 mm wide). Area d is an irregular pit (9.7 mm long x 6.8 mm 
wide) (Extended Data Fig. 1). 

The lower areas e-h are sets of multiple intersecting linear striations 
with V-shaped cross-sections, transverse to the long axis (maximum 
length 17.3 mm x 1 mm width) (Extended Data Fig. 1). Area i striations 
are sinuous, parallel (42.8 mm long) and transverse to the grain with 
indeterminate cross-sections. A small entry facet cuts area e. Area 
jis defined by multiple convex-shaping facets (27 mm maximum 
length x 14 mm width) intersecting the bark surface and forming the 
flat surface below (Extended Data Fig. 2). All facets have lost definition 
during storage. Area k (Extended Data Fig. 2) is a single convex facet 
(41.8 mm x 39.8 mm). 

On the left side and upper surface of part 2 (Fig. 4), at least 12 
convex-shaping facets and some intercutting are clearly visible on 
excavation (Extended Data Fig. 3). Five remain visible after storage as 
dimpling (Extended Data Fig. 3), resembling faceting marks on water- 
logged wood from Star Carr, UK”. Consecutive facets define a portion 
of the upper tapered edge of part 3 (Extended Data Fig. 4). Two clusters 
of small convex facets occur on the upper surface of part 1 (Fig. 4 and 
Extended Data Fig. 4). 

Surface modification occurs on the underlying trunk at its 
midpoint and along the narrowed end that passes through and beyond 
the notch (Fig. 3 and Extended Data Fig. 5). At the midpoint, a small 
area (approximately 10 cm x 5cm) preserves multiple short parallel 
striations transverse to the grain, with V-shaped sections, indicative of 
scraping’ (Extended Data Fig. 1). The log’s end tapers in plan and profile 
view with two sets of shaping marks visible. Parallel linear V-shaped 
cutmarks, transverse to the grain occur at the start of a breakin slope, 
increasing in depth and length downslope (Extended Data Fig. 5). 
Asecond set of marks occurs on the flat surface which emerges beyond 
the notch and continues into the section (Extended Data Fig. 5). Multiple 
linear and intersecting groups of V-shaped fine striations extend across 
this surface at acute and right angles to the grain. 

We interpret the notchas intentional, made by scraping and adzing' 
to create a join between the log and trunk, forming a construction of 
two connected parts. Infrared spectroscopy (Supplementary Infor- 
mation Section 5) provides indeterminate evidence for use of fire in 
shaping the notch. Clark” described a similar find, from the Acheulean 
in Site B, of comparable length (165 cm long) with a “wide and deepish 
groove’ transverse to the long axis, with tapered ends. He interpreted 
the groove as anthropogenic and suggested it was part of a structure. 


Fig. 3 | Structural unit formed by two overlapping logs in BLB5. The underlying log passes through a central notch cut into the upper log (object 1033) and 
extends into the section. Plan view of the unit (left) and during excavation (right). The numbers refer to the distance in centimetres. 
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Fig. 4 | Annotated images of the BLBS5 upper log (object 1033) showing areas 
of intentional modification. From left to right, the location of the central 
notch in profile, shaping marks in and on the margins of the notch (a-k), the 
notchin profile from the opposite side. The image on the right shows the upper 


The excavation of two interlocking logs in BLBS, with shaped ends on 
both objects supports this interpretation. We know of no comparable 
construction in the early archaeological record. 

In BLB3, two wood objects were found among plant material beside 
alarge unmodified treetrunk. Flake tools, cleavers, handaxes and core 
axes occur in this deposit. Object 661 (Ficus spp.) isa V-shaped piece of 
branch, lenticular in cross-section with no visible shaping marks (not 
illustrated). Object 660 (Kigelia africana), 36.2 cm long, rounded on 
one side with outer bark on both surfaces, tapers to an offset point cut 
60° across the long axis (Fig. 2b and Extended Data Fig. 6). The base is 
convex witha deep crack initiated in the pith, extending partially down 
each side (Extended Data Fig 6). Two modified areas occur below and 
transverse to the grain: a facet entry heel” (a) and faint parallel linear 
striations of indeterminate profile (b) (Extended Data Fig. 6). Above 
the tip are four sets of marks: two intercutting convex facets with inter- 
nal linear V-shaped cross-sections (c and d) and two sets of parallel 
striations (V-shaped) transverse to the grain (e and f) (Extended Data 
Fig. 6). The tip ends ina blunt break, subsequently rounded, its other 
side having three faint step fractures transverse to the grain obscured 
by the pitted surface (Extended Data Fig. 6). 

We interpret this as an intentionally shaped object altered by 
high-impact compression down its long axis. Possible uses include a 
wedge, an incompletely processed firewood and a portable work base 
(Methods). Minimally modified bark wedges are known ethnographi- 
cally’ and short-shaped roundwood wedges from the Mesolithic (Star 
Carr, UK)”. This unusual object has no parallel in the Middle Pleisto- 
cene record. 

In BLB2, object 219 (K. africana) was recovered in two refitting 
parts, one preserved in sediment, the other in the river directly below 
(Extended Data Fig. 7). The combined length is 62.4 cm, amaximum 
breadth at the base 11.9 cm, midsection breadth 6.1 cm and 1.3 cm at 
the offset tip (Fig. 2c). No bark occurs on the object; the upper surface 
is concave from base to tip (Extended Data Fig. 7). Faint linear stria- 
tions, with indeterminate cross-sections, cut across the grain near the 
tip on both sides along with small (less than 10 mm) indistinct convex 
facets (Extended Data Fig. 7). The obverse surface preserves four areas 
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surface of the log, and the three parts of the log (1-3) separated by cracks. 
White arrows indicate locations of shaping facets on the sides and upper 
surface of the log. 


with linear striations above the tip, including around two knots, with 
indeterminate cross-sections (Extended Data Fig. 7). Tip rounding 
may be from exposure to the river or use. Mid-Holocene waterlogged 
deposits at Gwisho, central Zambia preserved four digging sticks. Only 
one has small facets (less than 20 mm) on both sides of the tip”. The 
Gwisho tools are short by comparison with those used historically by 
Kalahari foragers” but similar in length to object 219. We interpret 
object 219 as a digging stick based on its morphology and presence 
of shaping facets®*. 

In BLB4 two wood objects were preserved in clays, one above 
the other (Fig. 2d,e). No stone tools were found. The upper object 
(Combretum zeyheri) is rectangular (59.24 cm x 29.34 cm x 7.7 cm) with 
bark traces and sapwood exposed on longitudinal surfaces (Extended 
Data Fig. 8). Evidence of radial cellular compression indicates flattening 
by sediment overburden. The ends bear distinct chop marks across the 
long axis (Fig. 2d; Extended Data Fig. 8). The obverse surface has one 
entry facet parallel to the grain and a set of transverse parallel linear 
striations of indeterminate profile (Extended Data Fig. 8). 

Well-preserved chop marks occur at enda where they descend in four 
steps, with a V-shaped profile, maximum depth of 2.3 cm, and the top 
step expanding in width from 13 mm to 25 mm (Extended Data Fig. 8). 
Stepped cuts occur across end b (Extended Data Fig. 8). The length of 
cuts and slight curvature suggest a broad, sharp, cleaver-like edge, 
the depth indicating applied direct force, either handheld or hafted 
(Methods and Supplementary Information Section 7). We interpret 
this object as a portion of treetrunk cut to size, indicating capacity to 
work wood at a large scale. 

The lower object (C. zeyheri) (Fig. 2e) is from branch wood, witha 
side-branch attachment, split centrally, length 37.9 cm, tapering from 
the base (12.3 cm wide) to a broken tip (2.1 cm wide) (Extended Data 
Fig. 9 and Supplementary Information Section 4). A single chop mark 
above the tip is transverse to the grain, V-shaped in profile with an 
irregular termination. We were unable to interpret this object. 

The recovery of modified wood with an exceptional level of preserva- 
tion at Kalambo Falls places woodworking firmly in the Mid-Pleistocene 
of Africa. Our understanding of innovation among early hominins has 


been largely inferred from lithic artefacts”. Wood from tree trunks 
enabled humans to construct large objects such as platform founda- 
tions that necessitated tools for felling and hewing. The large cutting 
tools of the Acheulean fulfilled these roles in the long co-evolution of 
stone- and woodworking” ©. Hominin use of the Kalambo River basin 
coincides with an extended period of forest cover (470-274 kyr)**””. 
Availability of forest resources and a permanent elevated water table 
created a habitat conducive to sustained occupation. Life in a peri- 
odically wet floodplain would be enhanced by constructing a raised 
platform, walkway or foundation for dwellings®. 

Woodworking was also a precursor for the Mid-Pleistocene invention 
of hafting, which increased efficiency of basic actions suchas chopping 
and scraping”. The interlocking logs from BLBS anticipate hafting’s 
core concept: the combination of two or more parts to make a con- 
struction”, enhancing our understanding of the technical cognition of 
these toolmakers”. Exceptional conditions of preservation give us this 
glimpse ofa capacity to create a built environment by hominins hitherto 
perceived as mobile foragers with limited technological diversity*. 
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Methods 


Sedimentology and depositional environments 

Sediment details were recorded from BLB1, BLB3 and BLBS (Fig. 1b) 
which focused on the oldest sediments exposed at the base of the river 
cliff. BLB3 and BLBS sampled the lowest 1-2 m of the cliff exposure; 
BLB1 spanned the full exposure from the base to the cliff top. A single 
sediment log of the BLB1 succession is shown (Fig. 1c) with only the 
lowest section reported in detail. 

Sediments are dominated by well-sorted sands exhibiting weakly 
developed cross-bedding up to 50 cm thick. Sands vary between fine/ 
medium to coarse grained. Gravels are rare, occurring as small discrete 
clasts in sandy units, or as lags associated with erosional horizons. Sedi- 
ments are consistent with lateral and downstream accretion of point 
and mid-channel bars within a sand-dominated fluvial channel system. 
Fine-grained sediments are less abundant but occur as either thin (less 
than 10 cm) discrete lenses or organic-rich clays and silts, frequently 
associated with fossil tree trunks or more laterally continuous units of 
silts and clays, their geometry indicative of accumulation in abandoned 
channels. The occurrence and arrangement of these differing units 
indicate a sandbed river that underwent regular low-flow conditions 
and episodically shifted, or avulsed, across its floodplain. Nothing in 
these sequences indicates major change in river process or substantial 
environmental transitions. Particle sizes were analysed using a Malvern 
Master Sizer 3000 (Supplementary Information Section 1). 


Luminescence dating 

Sixteen samples for dating were collected at Site BLB by hammering 
opaque plastic tubes into the sediment. A combination of field gamma 
spectrometry, laboratory alpha and beta counting and geochemi- 
cal analyses were used to determine radionuclide content, and the 
dose rate and age calculator” was used to calculate radiation dose 
rate. Sand-sized grains (approximately 150 to 250 um in diameter) 
of quartz and potassium-rich feldspar were isolated under red-light 
conditions for luminescence measurements and measured on Ris@ TL/ 
OSL instruments* using single-aliquot regenerative dose protocols***. 
Single-grain quartz OSL measurements dated sediments younger 
than around 60 kyr, but beyond this age the OSL signal was saturated. 
pIRIRSL measurements of aliquots consisting of around 50 grains of 
potassium-rich feldspars were able to provide ages for all samples 
collected. The pIRIRSL signal yielded an average value for anomalous 
fading of 1.46 + 0.50% per decade. Where quartz OSL and feldspar pIR 
IRSL were applied to the same samples, the ages were consistent within 
uncertainties without needing to correct for anomalous fading. The 
conservative approach taken here has been to use ages without any 
correction for fading. If a fading correction had been applied then the 
ages for the wooden artefacts would be older (Supplementary Informa- 
tion Section 2 and Supplementary Tables 2.8 and 2.9). 

Using the pIRIRSL ages without any fading correction, values at each 
of the five sections (BLB1 to BLB5; Extended Data Table 1) are in strati- 
graphic order (considering the calculated lo uncertainties on the ages), 
except for two samples. At BLB2, the central ages for samples KF10 
(355 + 18 kyr) and KF11 (422 + 28 kyr) are reversed. The reason for this is 
unknown. The age of this section is estimated by taking the mean of the 
two ages, giving a value of 383 + 28 kyr. The ages fall into four groups, 
containing between three and six samples in each group. The oldest 
cluster of ages contains three samples found at or below current river 
level in sections BLB3 (KF14: 492 + 32 kyr; KF15: 453 + 27 kyr) and BLBS 
(KFO1: 487 + 28 kyr). The second cluster is stratigraphically above the 
first and contains three samples, one in section BLBS (KFO2: 416 + 37 kyr) 
and the two ages at BLB2 that are in reverse order (KF10: 355 + 18 kyr 
and KF11: 422 + 28 kyr). The third cluster is stratigraphically above the 
second and contains four samples in BLB1 (KFO9: 355 + 33 kyr), BLB3 
(KF13: 326 + 16 kyr) and BLB4 (KF16: 322 + 19 kyr; KF17:316 + 18 kyr). The 
youngest group contains six samples found in BLB1 (KFO4: 11 +1 kyr; 


KFOS: 15 + 2 kyr; KFO6: 29 + 2 kyr; KFO7: 38 + 3 kyr; KFO8: 40 + 2 kyr) 
and BLB3 (KF12: 62 + 4 kyr). The ages of the three oldest clusters were 
calculated by combining ages using the central age model“, using only 
the randomerrors for weighting (Supplementary Information Section 2 
and Supplementary Tables 2.8 and 2.9), and adding in quadrature the 
systematic errors to determine the total uncertainty on the weighted 
mean ages. The ages for the clusters are 476 + 23 kyr, 390 + 25 kyr and 
324 +15 kyr. Ifa correction for anomalous fading of 1.46 + 0.50% per 
decade is applied, the ages for these three clusters and the associated 
artefacts are older (554 + 32 kyr, 452 + 29 kyr and 377 + 21 kyr). 


Wood identification and anatomy 

Subsamples of the specimens excavated from BLB2, BLB3, BLB4 and 
BLBS were examined in the biological preparation and scanning elec- 
tron microscope (SEM) laboratories in the Department of Scientific 
Research, British Museum. Given the three-dimensional (3D) nature 
of wood anatomy, each wood sample, irrespective of its size, was frac- 
tured manually to showtransverse, radial longitudinal and tangential 
longitudinal sections”. Each uncoated transverse, radial longitudinal 
and tangential longitudinal wood section was mounted onto SEM alu- 
minium stubs. 

SEM examination of the wood samples and comparative reference 
specimens (prepared and mounted using the same method) was under- 
taken ina variable pressure SEM, Hitachi S-3700N using the backscatter 
electron detector at 15 kV, with the SEM chamber partially evacuated 
(40 Pa). Magnifications ranged from x50 to x1,000. Preferred working 
distance was approximately 14 mm but was raised or lowered from 
12 mm to 16 mm as required. With the backscatter electron detector, 
3D mode (rather than compositional) was selected for maximum 
topographical information and to reveal diagnostic features for 
identification to genus or species level (Supplementary Information 
Sections 3 and 4). 


Wood surface analysis 
The BLB wood was airfreighted from Zambia to the University of Liv- 
erpool for analysis (2019-2022). Specialist advice was followed and 
the samples kept submerged in cold tap water in the dark with trace 
chlorine providing protection from fungal contamination. The water 
was changed regularly and the samples monitored for fungal growth; 
none was found. Definition of surface shaping marks on pieces from 
BLB3 (object 660) and BLBS (object 1033) degraded during storage, 
but no changes were seen on wood from BLB2 (object 219) and BLB4. 
The wood from the latter two contexts was recovered from clay and 
silt deposits with a high organic content, whereas the older contexts 
of BLB3 and BLBS were primarily sands. The lower piece from BLB4 
was found cracked in situ. Object 1033, BLBS, had two cracks (Fig. 2a) 
on either side of the central notch. The lower trunk was left in situ; 
excavation pits were backfilled. 

Samples are being conserved at York Archaeological Trust, UK, before 
being returned to Zambia for curation (Livingstone Museum). 

Radiocarbon dating and species identification samples were col- 
lected using a Chroma Gesellschaft F3400 Cork borer (5 mm internal 
diameter). Sample depth ranged from 8 mm to 12 mm. Samples for 
infrared and Raman spectroscopy (to identify mineralization and 
burning) were collected from BLBS (object 1033) and BLB2 (object 
219) using a Swann-Morton no. 3 scalpel with 10a blades. Samples 
were 5mm x 5mm x 2 mm or less and freeze dried before analysis to 
remove moisture. Spectra were recorded using samples larger than 
attenuated total reflectance diamond crystals. Opus 7.2 software was 
used to process sample spectra using 30 scans covering a wavelength 
range of 500-4,000 cm”. An infrared spectrum of modern C. zeyheri 
was provided by Royal Botanic Garden Kew as a comparator, with 
spectra recorded on three separate regions. Silica peaks occurred in 
the archaeological samples but were absent in the modern samples, 
indicating mineralization of the BLB wood. Asample collected fromthe 


notch in BLBS5 (object 1033) showed a flattened silica peak suggestive 
of carbonization from burning”. No evidence of burning was found 
using portable X-ray fluorescence (pXRF) and SEM energy-dispersive 
X-ray spectroscopy analyses of samples from BLBS (object 1033), BLB3 
(object 660) and BLB2 (object 219). Raman spectroscopy was unsuc- 
cessful; samples were obliterated by the laser before analyses were 
completed (Supplementary Information Section 5). 

Allwood was photographically modelled while submerged to mini- 
mize damage from drying (University of Liverpool). RTI was used to 
produce two-dimensional models/images of complete specimens. 
RTI models differ from a standard photographic image in being view- 
able from different angles, allowing all lighting angles to be observed 
as well as exposing low-relief surface features. RTI photographs were 
taken using a Sony Alpha 77ii DSLR using manual settings—F stop value 
F16, witha shutter speed of 1/200—ISO-200. During photography, the 
wood was submerged and placed close to the surface of the water to 
minimize distortion. RTI models were built using RTI builder v.2.0.2, 
HSH highlight-based fitter v.1.0.1, and viewed with RTI viewer v.1.1.0. 
(www.Culturalheritageimaging.org). Measurements from RTI models 
were taken with image (v.1.53t) and TPsDig2 v.2.31. GIMP v.2.10 was 
used for postprocessing RTI images (removal of background and inser- 
tion of digital scales). Underwater photogrammetry was performed 
using acustom-built watertight container housing acamera. Two cam- 
eras were used depending on the size of object being modelled: Sony 
A77M2 with a Sigma 30 mm lens for larger objects and Sony A7M3 
with a Sony 90 mm macro lens for smaller objects. Photographs were 
taken in raw format and converted to JPEG using Adobe Lightroom. 
Three-dimensional photogrammetry models were constructed using 
Agisoft Metashape (v.1.8.1, build 13915, 64 bit) selecting the highest 
accuracy settings. All 3D models were refitted and rotated for imaging 
as.stl files using MorphoDig (v.1.6.7 64 bit). Additional photographs 
were taken after objects were allowed to dry for upto 15 min to increase 
visibility of surface features. 


Water modification and saturation experiments 

A pilot study was undertaken to assess the effects of flowing water 
and water saturation on preservation of modified wood surfaces. The 
preservation of shaping marks on wood at Site BLB contrasts with the 
reported absence of surface marks on wood recovered 60 m down- 
stream at Site B, in similar late Acheulean deposits”. Site B is on the 
actively eroding outer bend of the Kalambo River, whereas BLB lies on 
astraight stretch of the main channel (Fig. 1b). Controlled experimen- 
tation was undertaken involving modification of standardized wood 
samples using stone tools and a steel knife, with the surfaces recorded 
before and after a period of exposure to abrasive flows. Thirty-two 
specimens of soft wood (Abies alba, silver fir) of standardized size were 
hand-stripped of their outer bark, eight specimens were left unmodified 
as a control sample, four were allowed to dry before placing in the flume 
tank and four were kept wet. The specimens analysed were branches 
taken from trees of similar ages and growth. A soft wood was chosen 
for its low-density fibre structure which increases the likelihood of 
modification by an abrasive flow regime and swelling from saturation. 
The 24 worked specimens were split into three groups of eight: one 
group was worked witha steel knife, one group with a bifacially worked 
flint tool and one using a unifacially retouched flake tool. A steel knife 
was incorporated to increase the contrast in edge morphologies and 
marks produced. All specimens were cut to similar sizes, with lengths 
ranging from 85.6 mm to 118 mm, and were scraped unidirectionally 
for 10-15 min until surfaces were flat to enhance the identification 
of minor modifications caused by abrasive flow and those made by 
intentional shaping. Analysis was primarily visual, including inspection 
with the naked eye, SEM analysis and RTI. The combination of analyti- 
cal techniques provided functional and morphological information. 
Blind tests were undertaken to assess the analyst’s ability to identify 
working methods and tools using RTI and SEM. 


A GUNT HM166 flume tank was used which gives a flow velocity of up 
to1.5ms. The speed was set at approximately 1 m s” (0 to 2.0 Knots) 
with water depth of 80 mm to cover the sticks and a channel width 
of 50 mm, resulting in a proportional flow rate equivalent to an aver- 
age river flow*’. The experimental channel was filled with coarse- and 
fine-grained sand, like those from the basal levels of Site BLB. The wood 
specimens were secured in a wooden rig and periodically rotated to 
ensure equal exposure of surface to the path of the flowing sediments. 
The flow velocity allowed for the greatest possible effect on the sticks, 
which were left in the water for 2 weeks, 70 hours of which involved 
exposure to the flow, and the remainder in still water overnight and 
weekends. In addition, four of each group of eight sticks were kept in 
water after exposure to the flume tank and before analysis to assess 
the effect that water absorption had on surface modification visibility. 
Once the 2-week period was completed, half the sticks were allowed to 
dry and were re-analysed using the methods above. The results from 
each set of analyses were compared, to assess the differences before 
and after exposure to flowing water and the effects of saturation on 
the preservation of surface morphology. 

The experiments showed that (1) abrasive flows smoothed facet 
edges and ridges, but tool marks remained visible; (2) saturated speci- 
mens lost definition of facet junctions from swelling, and minor surface 
irregularities from shaping were smoothed or disappeared altogether; 
and (3) fluvial exposure did not create marks mimicking intentional 
shaping marks. Differences between the four dry control samples and 
the four saturated samples were only visible under SEM examination 
(Supplementary Information Section 6). 


Woodworking experiments 

To characterize the patterns of wear left on wood from working and 
shaping with stone tools at Kalambo Falls, four transverse scraping 
experiments were performed at the University of Liverpool anda 
chopping experiment undertaken at the University of Liége, Traceo- 
Lab. Quartzite was used because of its prevalence as a raw material of 
choice for making small and large tools”. The quartzite was sourced 
inthe UK and from Kalambo Falls. Modern reference samples of wood 
found at Kalambo could not be sourced. Two European woods of 
similar dry density were selected for their comparability to the two 
species that feature in the Kalambo sample, C. zeyheri (0.58 g cm)” 
(BLB4, BLBS5) and K. Africana (0.58-0.64 g cm’’) (BLB2)*°. The Euro- 
pean species were ash (Fraxinus excelsior) (0.56-0.64 g cm) and elm 
(Ulmus procera and Ulmus minor) (0.55-0.60 g cm’*) (https://www. 
engineeringtoolbox.com/wood-density-d_40.html). Elm has long, 
intertwined fibres that give strength, making it a less tractable hard- 
wood for working with stone. We do not have comparable fibre length 
data for the Zambia woods. Retouched and un-retouched flakes were 
used unidirectionally, handheld, following the grain for 20 min each 
to scrape bark-covered green ash and green elm. A quartzite biface was 
hafted (juxtaposed)™ to a wood handle made by a traditional craftsman 
from Chiungu Village, Kalambo Falls and used in a scraping motion 
for 20 min. 

A single wedge experiment was undertaken using a tapered cut 
section of hazel to laterally split open a larger section of birch, result- 
ing in perpendicular compression damage on the active face of the 
wedge in contact with the birch. The wedge tip, which was initially not 
perpendicular, snapped during use leaving a perpendicular break. 
The struck end of the wedge displays a series of small radial cracks. 
The asymmetric and undulating morphology of the tapered end of 
BLB3 (object 660) (Fig. 2b and Extended Data Fig. 6) was replicated 
in experimental burning and scraping of hazel wedge tips using flint 
flakes. No other evidence exists for burning on object 660. The deep 
crack on the base of object 660 (Extended Data Fig. 6) was replicated 
by a splitting experiment using a flint handaxe placed point down on 
the butt of a hazel pole and struck using a wood baton. A crack propa- 
gated from the point impact and spread one-third down the pole. The 
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cracking of object 660 may have resulted from intentional splitting to 
create firewood, or from use as a portable work surface (‘stake anvil’), 
with a stone working edge inserted into the base enabling bimanual 
processing of organic materials (stake set into firm ground and gripped 
between knees for support). 

The chop marks visible at end a of BLB4 upper object (compressed 
log) (Extended Data Fig. 8) were replicated by chopping fresh elm 
using a quartzite cleaver hafted in a socketed wooden handle (elm). 
The tool proved effective and succeeded in producing the deep cuts 
seen ona, and the distinctive steps on a were produced automatically 
during chopping. Depth of the cuts varies between 5 mm and 8 mm, 
and their width corresponds to the width of cleaver used (77 mm). The 
cleaver remained sharp and effective for about 15-20 min of work, after 
which the edge became blunt and tore fibres more than it was cutting 
them. A blunter tool edge also results ina more fibrous appearance of 
the chopped wood. The depth of cuts reproduced experimentally is 
slightly shallower than the 13-25 mm visible on the BLB4 log. However, 
the tool used was lighter than archaeological examples found at the 
site (experimental cleaver: 402 g, with handle 1,020 g; which is well 
below the 1,787 g of an archaeological example from BLBS), its spine 
plane angle more obtuse (30-35° compared with the 29-31° of the 
archaeological cleaver) and its edge already showed some removals 
before use. Cuts of 13-15 mm in depth, as on the BLB4 log, are eas- 
ily attainable with a heavier cleaver with a more acute edge angle 
(Supplementary Information Section 7). 


Ethnographic observations in Zambia 

A traditional woodworker, John Mukopa, Mungwi area (Kasama), 
Northern Province, participated in ethnographic research undertaken 
26-28 July 2022 by Moto Moto Museum, and University of Liverpool 
staff were supported by the Endangered Material Knowledge Pro- 
gramme (https://www.emkp.org) (EMKP2020SG01). Ethics approval 
was obtained from the University of Zambia (HSSREC-20220-APR-022), 
the University of Liverpool (REF: 11408) and the British Museum/EMKP. 
Observations were made onthe trees selected, tools used in felling and 
shaping, and marks left on wood surfaces from shaping. Inferences 
were drawn for interpreting the production sequence in making object 
1033 (BLBS) and possible uses of the log with chopped ends (BLB4) 
(Supplementary Information Section 8). 


Reporting summary 
Further information on research design is available in the Nature Port- 
folio Reporting Summary linked to this article. 


Data availability 


Photographic data (wood objects, RTI and photogrammetry models) 
supporting this research are openly available via Liverpool Research 
Data Catalogue (https://doi.org/10.17638/datacat.liverpool.ac.uk/ 
2200). 
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Extended Data Fig. 1| Location of modification marks on object 1033 (BLB5). and impact pit ‘d’ (scale =3 cm). Lower: 2a shows striation clusters ‘e’ and f’ 
Twoareas of modification are identified in the Upper and Lower parts of the (scale =3 cm). 2b. Intersecting striations at ‘g’; parallel striations at ‘h’; wavy 
notch (Fig. 2a; Fig. 4). Upper: 1a shows distribution of marks ‘a-d’; 1b provides bands of flat striations at ‘i’ (scale=3 cm). 

detail of intersecting striations ‘a’ and facet ‘b’. Image 1c shows facets ‘b’ and ‘c’ 
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Extended Data Fig. 2| Location of modification marks on Part 1, object highlights the loss of facet definition in area ‘j’ during storage. (Scale =3 cm). 
1033. Upper left image shows area ‘j’ on the lower part of the notch seen on The fine surface detail is gone leaving only the deepest scarring. Lower image 
excavation with distinctive and extensive convex chop marks onthe surface is‘k’,asingle convex chop-mark on lower side of the notch, cutting into the 


and back edge where it intersects the upper surface of the log. Right image surface bark. (Scale =1 cm). 


Extended Data Fig. 3 | Shaping marks on Part 2 of object 1033. Upper image 
shows shaping marks on the tapered side and surface of Part 2 as seen on 
excavation (Fig. 4). Lower image shows loss of definition during storage. The 
arrows in the Lower image correspond with the area to the left of centre onthe 
Upperimage. 
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Extended Data Fig. 4 | Shaping marks on the upper surfaces of object 1033 
and onthe underlying treetrunk (Fig. 3). Clockwise, from left upper left; 
chop marks on Part 2 (Extended Data Fig. 3); cluster of small convex hewing 
marks on Part 1, near Part 2 (Fig. 4); cutmark (upper arrow) and small facets 
(lower arrows) on Part 1 near Part 3 (Fig. 4); intercutting chop marks onthe 


upper right edge of the Part 3 taper (Fig. 4); underlying log midsection, 
intersecting cutmarks transverse to the grain (bold arrow, upper left, 
indicating direction of grain). Marks on underlying treetrunk interpreted 
as result of scraping, perhaps from debarking. 


Extended Data Fig. 5 | Shaping marks on the treetrunk underlying object 
1033 notch, BLBS. Surface modifications on the treetrunk (Fig. 3) beneath and 
just beyond the notch (Fig. 3., Extended Data Fig. 4). Images described counter 
clockwise: white arrows highlight direction of striations, some V-shapedin 
profile, cutting across the grain from the edge toward the middle of the 
surface; aset of parallel transverse striations/cutmarks, with V-shaped profiles 


(white arrows) along ridge descending toa flat tapered surface; inset detail 

of the striations on the ridge crest; overview of the set of striations along the 
ridge; profile view of the flattened surface below the ridge; profile view of the 
ridge and break of slope onto the flat tapering surface which extends beyond 
the notch; plan view of the flat tapering end with transverse striations and 
cutmarks with V-shaped profiles (white arrows) interpreted as shaping marks. 
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Extended Data Fig. 6 | Shaping and possible use marks on object 660, BLB3. damage, and ‘b’ indicates a set of faint transverse striations; ‘c’ and ‘d’ 
Surface modifications on ‘wedge’ (Fig. 2b) shown clockwise: white arrows intersecting facets above tip, ‘e’ and ‘f’ are sets of parallel faint transvers 
indicating location of modifications on base and either side of the tip; basal striations; tip face with convex step terminations, possibly from use. 
crack (top arrow) and side split, facet ‘a’ intercepts a striation, possible use 


Extended Data Fig. 7 | Shaping marks on object 219, BLB2 (Fig. 2c). Upper facets away from the tip. (Scale =1 cm.) The vertical image to the right shows 
image shows all views of the object with white arrows indicating location of four groups of faint parallel transverse striations in an area with two knots, 
marks oneither side of the tip and above the tip. Both sides of the tip (two noscale. The object is interpreted as a digging-stick. 

images on left) havea faint set of transverse striations and two or more small 


Article 


Extended Data Fig. 8| Chop marks on ‘cut log’, BLB4 (Fig. 2d). Upper image width of the step between the first and second chop-mark expanding, left to 
shows all views of the object and location of modifications (white arrows) on right, from 13 to 25 mm. The deep, slightly convex cut of the first chop is visible. 
the ends, ‘a’ and ‘b’, two marks on the underlying surface and the preservation The underly step expands left to right from 6 mmto1l mm. A third step starts 
of bark on both surfaces. Image below is ‘a’ with its distinctive stepped chop witha width of 7 mm. Bottom image is end ‘b’ with black arrows indicating 


marks (black arrows), and underlying grey image (denatured RTI) shows the location of chop marks. 


Extended Data Fig. 9 | Images of split branch, BLB4 (‘notched stick’; Fig. 2e). Above, close-up of the facet/chop-mark with arrow indicating direction of 
initiation. Below, all views of split branch with chop (white arrow) and close-up of the entry facet. The object was found cracked in situ. 
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Extended Data Table 1| Luminescence ages at Site BLB, Kalambo Falls. Six samples could be dated using the quartz optically 
stimulated luminescence signal (OSL). All 16 samples were dated using the post-infrared infrared stimulated luminescence 
(post-IR IRSL) signal from potassium-rich feldspar grains 


Quartz OSL Feldspar Post-IR IRSL 
Site Sample De* Dose Rate Age t Det Dose Rate Age t 
(Aber232/) (Gy) (Gy/kyr) (kyr) (Gy) (Gy/kyr) (kyr) 
BLBS KFO1 859 + 29 1.76 + 0.08 487 + 28 
KFO2 884 + 49 2.13 +0.15 416 + 37 
BLB1 KF04 1841 1.65 + 0.07 1141 29+1 2.59 + 0.11 1141 
KFO5 1841 1.51 + 0.07 1241 36 +3 2.36 + 0.16 1542 
KFO6 38 +1 1.32 + 0.06 28 +2 63 +2 2.15 +0.15 2942 
KFO7 3742 0.97 + 0.04 3742 69+2 1.79 + 0.15 38 +3 
KF08 552 1.51 + 0.07 36 +2 94 +4 2.36 + 0.10 4042 
KFO9 634 + 29 1.79 + 0.15 355 + 33 
BLB2 KF10 676 + 18 1.91 + 0.09 355 + 18 
KF11 818 + 39 1.94 + 0.09 422 + 28 
BLB3 KF12 56 +2 1.34 + 0.06 4143 13447 2.15 + 0.09 62+4 
KF13 587 +11 1.80 + 0.08 326 + 16 
KF14 783 + 36 1.59 + 0.07 492 + 32 
KF15 758 + 26 1.67 + 0.08 453 + 27 
BLB4 KF16 714 +10 2.22 + 0.13 322 +19 
KF17 724 +19 2.29+0.11 316 + 18 


* Six hundred quartz grains of each sample had their OSL signal measured, and between 86 and 176 yielded equivalent dose (De) values. These were combined 
using the central age model (CAM*°) 


t Uncertainties shown for the ages listed here (shown at 1 — sigma) include both random and systematic sources of uncertainty. 


+ Twelve aliquots of potassium-rich feldspar were measured to determine De for each sample, and all passed the acceptance criteria except one aliquot of KF 12. 
The CAM was used to calculate the weighted mean dose. 
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Statistics 
For all statistical analyses, confirm that the following items are present in the figure legend, table legend, main text, or Methods section. 


n/a | Confirmed 


| The exact sample size (n) for each experimental group/condition, given as a discrete number and unit of measurement 


| A statement on whether measurements were taken from distinct samples or whether the same sample was measured repeatedly 


The statistical test(s) used AND whether they are one- or two-sided 
Only common tests should be described solely by name; describe more complex techniques in the Methods section. 


A description of all covariates tested 


| A description of any assumptions or corrections, such as tests of normality and adjustment for multiple comparisons 


E A full description of the statistical parameters including central tendency (e.g. means) or other basic estimates (e.g. regression coefficient) 
AND variation (e.g. standard deviation) or associated estimates of uncertainty (e.g. confidence intervals) 


m For null hypothesis testing, the test statistic (e.g. F, t, r) with confidence intervals, effect sizes, degrees of freedom and P value noted 
Give P values as exact values whenever suitable. 


For Bayesian analysis, information on the choice of priors and Markov chain Monte Carlo settings 
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Estimates of effect sizes (e.g. Cohen's d, Pearson's r), indicating how they were calculated 


Our web collection on statistics for biologists contains articles on many of the points above. 


Software and code 


Policy information about availability of computer code 


Data collection Reflectance transformation imaging (RTI) models of the wood objects were built using RTIbuilder 2.0.2, HSH highlight based fitter 1.0.1, and 
viewed with RTI viewer 1.1.0. (www.Culturalheritageimaging.org). Measurements from RTI models were taken with imageJ (version 1.53t) 
TPsDig2 version 2.31. GIMP 2.10 was used for post processing RTI images. 
Photographs of the wood were taken in raw format and converted to JPEG using Adobe Lightroom. 3D photogrammetry models were 
constructed using Agisoft Metashape (1.8.1, build 13915, 64 bit) selecting the highest accuracy settings. 3D models were refitted and rotated 
for imaging as .stl files using MorphoDig (version 1.6.7 64 bit). 


Data analysis Opus 7.2 software was used to process the infrared spectra of each waterlogged wood sample. 


For manuscripts utilizing custom algorithms or software that are central to the research but not yet described in published literature, software must be made available to editors and 
reviewers. We strongly encourage code deposition in a community repository (e.g. GitHub). See the Nature Portfolio guidelines for submitting code & software for further information. 


=> 
jad) 
v 
= 
= 
D 
© 
O 
Be 
2 
= 
= 
D 
© 
fe) 
=; 
D 
a 
Za) 
= 
3 
3 
w% 
= 
= 


Data 


Policy information about availability of data 
All manuscripts must include a data availability statement. This statement should provide the following information, where applicable: 


- Accession codes, unique identifiers, or web links for publicly available datasets 
- A description of any restrictions on data availability 


- For clinical datasets or third party data, please ensure that the statement adheres to our policy 


Photographic data (wood objects, RTI models, and 3D photogrammetry models) freely available (open access) at 10.17638/datacat.liverpool.ac.uk/2200 


Research involving human participants, their data, or biological material 


Policy information about studies with human participants or human data. See also policy information about sex, gender (identity/presentation), 
and sexual orientation and race, ethnicity and racism. 
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Reporting on sex and gender Ethnographic research (Supplementary Section 8) was undertaken in Zambia with a single participant (self-reported as 
biological male, aged 80) with full informed consent given. 


Reporting on race, ethnicity, or Ethnicity self-reported as Bemba (linguistic identification). 
other socially relevant 


groupings 

Population characteristics The study is based on a single individual. 

Recruitment The individual was selected as one of the last three remaining traditional wood-working craftspeople in northern Zambia. 
Informed consent was obtained and ethics approval was obtained from the University of Zambia (HSSREC-20220-APR-022) 
and the University of Liverpool Central Univeristy Ethics Comitteees (REF: 11408). 

Ethics oversight Ethics approval was obtained from the University of Zambia (HSSREC-20220-APR-022), the University of Liverpool Central 
University Ethics Committees (REF: 11408) and the British Museum (EMKP 2020SG01) with end of project reporting to these 
institutions. 


Note that full information on the approval of the study protocol must also be provided in the manuscript. 


Field-specific reporting 


Please select the one below that is the best fit for your research. If you are not sure, read the appropriate sections before making your selection. 


[| Life sciences Behavioural & social sciences Ecological, evolutionary & environmental sciences 


For a reference copy of the document with all sections, see nature.com/documents/nr-reporting-summary-flat.pdf 


Life sciences study design 


All studies must disclose on these points even when the disclosure is negative. 


Sample size Describe how sample size was determined, detailing any statistical methods used to predetermine sample size OR if no sample-size calculation 
was performed, describe how sample sizes were chosen and provide a rationale for why these sample sizes are sufficient. 


Data exclusions | Describe any data exclusions. If no data were state so OR if data were excluded, describe the exclusions and the 


rationale behind them, indicating whether exclus 


cluded from the analyses 


on criteria were pre-established. 


Replication If all attempts at replication were successful, confirm this 


ibe why. 


Randomization ted into experimental groups. If allocation was not random, describe how covariates 


were controlled OR if this is not relevant to your study, explain why. 


cribe whether the i 


tigators were blinded to group allocation during data collection and/or analysis. If blinding was not possible, 


Blinding 


ribe why OR explain why blinding was not relevant to your study. 


Behavioural & social sciences study design 


All studies must disclose on these points even when the disclosure is negative. 


Study description 


Research sample 


Sampling strategy 


Data collection 


Timing 


Data exclusions 


Non-participation 


Randomization 


Briefly describe the study type including whether 


data are quantitative, qualitative, or mixed-methods (e.g. qualitative cross-sectional, 
quantitative experimental, mixed-methods case st 


udy). 


State the research sample (e.g. Harvard university undergraduates, villagers in rural India) and provide relevant demographic 


information (e.g. age, sex) and indicate whether the sample is representative. Provide a rationale for the study sample chosen. For 
studies involving existing datasets, please describe the dataset and source 


Describe the sampling procedure (e.g. random, snowball, stratified, convenience). Describe the statistical methods that were used to 
predetermine sample size OR if no sample- alculation was performed, describe how sample 
zes are sufficient. For qualitative data, please indicate whether data saturation was considered, and 
what criteria were used to decide that no further sampling was needed. 


zes were chosen and provide a 


rationale for why these sample 


Provide details about the data collection procedure, including the instruments or devices used to record the data (e.g. pen and paper, 
computer, eye tracker, video or audio equipment) whether anyone was present besides the participant(s) and the researcher, and 


whether the researcher was blind to experimental condition and/or the study hypothesis during data collection. 


Indicate the start and stop dates of data collection. If there is a gap between collection periods, state the dates for each sample 
cohort. 


If no data were excluded from the analyses, state so OR if data were excluded, provide the exact number of exclusions and the 
rationale behind them, indicating whether exclusion criteria were pre-established. 


State how many participants dropped out/declined participation and the reason(s) given OR provide response rate OR state that no 
participants dropped out/declined participation. 


If participants were not allocated into experimental groups, state so OR describe how participants were allocated to groups, and if 
allocation was not random, describe how covariates were controlled. 


Ecological, evolutionary & environmental sciences study design 


All studies must disclose on these points even when the disclosure is negative. 


Study description 


Research sample 


Sampling strategy 


Data collection 


Timing and spatial scale 


Data exclusions 


Reproducibility 


Randomization 


Blinding 


Archaeological excavations took place at site BLB, Kalambo Falls, Zambia from 7 July to 4 August 2019. Five distinct areas were 
excavated, labeled BLB1, BLB2, BLB3, BLB4, and BLB5. BLB1 was a geological section through the full 9m depth of deposit. 


Five fossil wood objects were collected for analysis and 16 sediment samples for dating. 


The placement of excavation units was designed to sample a stretch of exposed river bank for sampling sediments and recovery of 
archaeological and paleoenvironmental material. The deepest excavations (BLB3, BLB5) reached 1 m below the water table at which 
point the excavations were halted. 


All objects and sediment samples were recorded in three-dimensional space using a total station by the project team under the 
direction of PI Prof Barham. Excavators wore nitrile powder free gloves to minimise surface contact and all artefacts were doubled 
bagged and separate sediment sample included for residue analyses. Artefact bags numbered and logged. Photographs for 
photogrammetry taken and sections drawn. 


Data collection at Kalambo Falls took place between 7 July and 4 August 2019 with wood objects recovered from four of the five 
excavation blocks. Ethnographic fieldwork took place in northern Zambia in July-August 2022. Analyses of the wood objects began in 
October 2019 at the University of Liverpool and was interrupted by Covid-19 from March 2020 until January 2021. Photographic 
analyses completed in October 2022. Infrared analyses of wood samples completed January 2023. 


No data were excluded. 


The data collection parameters for producing wood infrared spectra differed slightly between the analysis of the modern reference 
collection and the waterlogged wood and the modern reference sample. Royal Botanical Gardens, Kew (RBG) supplied a modern 
comparison sample of Combretum zeyheri with IR spectra from three separate areas on the wood. There was no preparation of the 
samples before analysis. Four scans were performed for each sample in the Mid-IR (4000-550 cm-1) and 32 in the Near-IR 
investigations (1000-4000 cm-1). Six samples were taken from each waterlogged wood specimens (object 219, 1033), and freeze 
dried before analysis. The spectra of each sample spectra of each sample was captured using 30 scans covering a wavelength range 
of 4000-500 cm-1. 


Infrared samples were collected within and beyond the target area of interest on object 1033 (BLB5) to identify patterning in the 
spectra. No randomisation was applied. 


Blinding was not necessary for the collection of infrared data. A modern reference sample of the wood species Combretum zeyheri 
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Blinding was used as the comparator to the archaeological wood identified as Combretum zeyheri. 


Did the study involve field work? Yes No 


Field work, collection and transport 


Field conditions Exposed riverbank sediments, including waterlogged sediments below river level. 


Location 082 35' 38.96"South; 312 14' 36.77", East, elevation 1390 amsl, river depth unknown 


Access & import/export Permission to excavate Site BLB, Kalambo Falls was granted by the National Heritage Conservation Commission as part of the Deep 
Roots Project (Permit dated 29 June 2017; Pl Lawrence Barham, University of Liverpool). The right to temporary export was granted 
by the National Heritage Conservation Commission on the agreement that the wood objects would be conserved before their return 
to the Livingstone Museum, Zambia (Permit NHCC/2019/01/LB1). Access to the site is by track from the town of Mbala, 35km to the 
south. A footpath leaves the track toward the river and follows the river bank. Wood samples were collected from waterlogged 
sediments and from organic clays above the river. The samples were stored in waterfilled tanks on site before being wrapped in cling 
film, to retain moisture, and embedded in 60 mm foam to minimise shock on travel, and placed in purpose made marine plywood 
boxes. The crated samples were transported to Lusaka from where they were airfreighted to the University of Liverpool for storage 
in water filled tanks during the period of analysis. The samples are now being conversed by the York Archaeological Trust. 


Disturbance All excavation areas were backfilled to minimise disturbance to the deposits. 


Reporting for specific materials, systems and methods 


We require information from authors about some types of materials, experimental systems and methods used in many studies. Here, indicate whether each material, 
system or method listed is relevant to your study. If you are not sure if a list item applies to your research, read the appropriate section before selecting a response. 


Materials & experimental systems Methods 
n/a | Involved in the study n/a | Involved in the study 
l Antibodies ChIP-seq 
l Eukaryotic cell lines Flow cytometry 
Palaeontology and archaeology MRI-based neuroimaging 


Animals and other organisms 


Dual use research of concern 


Plants 


l 
l Clinical data 
l 
l 


Antibodies 


Antibodies used Describe all antibodies used in the study; as applicable, provide supplier name, catalog number, clone name, and lot number. 


Validation Describe the validation of each primary antibody for the species and application, noting any validation statements on the 
manufacturer's website, relevant citations, antibody profiles in online databases, or data provided in the manuscript 


Eukaryotic cell lines 


Policy information about cell lines and Sex and Gender in Research 


Cell line source(s) State the source of each cell line used and the sex of all primary cell lines and cells derived from human participants or 
vertebrate models. 


Authentication Describe the authentication procedures for each cell line used OR declare that none of the cell lines used were authenticated. 


Mycoplasma contamination Confirm that all cell lines tested negative for mycoplasma co 
mycoplasma contamination OR declare that the cell lines 


tamination OR describe the results of the testing for 


were not tested for mycoplasma contamination. 


Commonly misidentified lines Name any commonly misidentified cell lines used in the study and provide a rationale for their use. 
(See ICLAC register) 
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Palaeontology and Archaeology 


Specimen provenance All specimens were collected from the site BLB, Kalambo Falls, Zambia (082 35' 38.96"South; 312 14' 36.77) with specimens identified 
by excavation area (BLB2, BLB3, BLB4, BLB5) and position recorded in x,y,z coordinates. Permission to excavate Site BLB was granted 
by the National Heritage Conservation Commission as part of the Deep Roots Project (Permit dated 29 June 2017; PI Lawrence 
Barham, University of Liverpool). The right to temporary export of the wood was granted by the National Heritage Conservation 
Commission on the agreement that the wood objects would be conserved before their return to the Livingstone Museum, Zambia 
(Permit NHCC/2019/01/LB1). 


Specimen deposition The specimens are currently being conserved by the York Archaeological Trust (York, UK) and on completion of conservation they will 
be deposited in the Livingstone Museum, Livingstone, Zambia. 


Dating methods Sixteen samples for luminescence dating were collected by hammering opaque plastic tubes into the sediment. The ends of the tubes 
were sealed in the field to avoid light exposure and to retain moisture, and then returned to the laboratory for analysis. Sediment 
from the inner part of the tubes was isolated and further processed for luminescence dating analyses. The samples were airfreighted 
from Zambia to Aberystwyth University, UK. All sample processing and analysis was undertaken at Aberystwyth Luminescence 
Research Laboratory in red light conditions. Prior to chemical treatments, the coarsest samples were wet-sieved to isolate grains 
ranging from 90-1000 um, whilst other samples were immediately subjected to chemical reagents. Carbonates were removed by 
treating the samples with 10% HCI for at least 3 days or until no continued reaction was observed. Initial reactions with HCI were mild 
or unobservable for all samples. After rinsing three times with distilled water, the sediment was put in to 20 vol H202 to dissolve 
organic matter before the samples were again rinsed three times and dried. Samples were then dry sieved using a stack of sieves 
ranging from 90 to 355 um. A combination of field gamma spectrometry, laboratory alpha and beta counting, and geochemical 
analyses were used to determine the radionuclide content and DRAC was used to calculate the radiation dose rate. Sand-sized grains 
(~150 to 250 um in diameter) of quartz and potassium-rich feldspar were isolated under red-light conditions for luminescence 
measurements and measured on Ris@ TL/OSL instruments using single aliquot regenerative dose protocols. Single grain quartz 
optically stimulated luminescence (OSL) measurements were able to date sediments younger than ~60 ka; beyond this age the OSL 
signal was saturated. Post-infrared infrared stimulated luminescence (pIR-IRSL) measurements of aliquots consisting of ~50 grains of 
potassium-rich feldspars were able to provide ages for all samples collected. The pIR-IRSL signal yielded an average value for 
anomalous fading of 1.46 + 0.50 % per decade. Where quartz OSL and feldspar pIR-IRSL were applied to the same samples, the ages 
were consistent within uncertainties without needing to correct for anomalous fading. 
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Tick this box to confirm that the raw and calibrated dates are available in the paper or in Supplementary Information. 


Ethics oversight Ethics oversight was required on the collection of ethnographic observations as outlined above, and the research was approved by 
the University of Liverpool Central University Ethics Committees (no. 11408) the University of Zambia Human and Social Sciences 
Research Ethics Committee (HSREC-2022-APR-022,and the British Museum Endangered Material Knowledge Programm (2020SG01) 


Note that full information on the approval of the study protocol must also be provided in the manuscript. 


Animals and other research organisms 


Policy information about studies involving animals; ARRIVE guidelines recommended for reporting animal research, and Sex and Gender in 
Research 


Laboratory animals For laboratory animals, report species, strain and age OR state that the study did not involve laboratory animals 


Wild animals Provide det 


Reporting on sex 


Field-collected samples 


Ethics oversight 


Note that full information on the approval of the study protocol must also be provided in the manuscript. 


Clinical data 


Policy information about clinical studies 
All manuscripts should comply with the ICMJE guidelines for publication of clinical research and a completed CONSORT checklist must be included with all submissions. 


Clinical trial registration | Provide the trial registration number from ClinicalTrials.gov or an equivalent agency. 


Study protocol Note where the full trial protocol can be accessed OR if not available, explain why. 
Data collection Describe the settings and locales of data collection, noting the time periods of recruitment and data collection. 


Outcomes Describe how you pre-defined primary and secondary outcome measures and how you assessed these measures. 


Dual use research of concern 


Policy information about dual use research of concern 


Hazards 


Could the accidental, deliberate or reckless misuse of agents or technologies generated in the work, or the application of information presented 
in the manuscript, pose a threat to: 


No | Yes 
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Public health 


|__| National security 


Crops and/or livestock 


Ecosystems 


Any other significant area 


Experiments of concern 


Does the work involve any of these experiments of concern: 


No | Yes 


Demonstrate how to render a vaccine ineffective 


Confer resistance to therapeutically useful antibiotics or antiviral agents 


Enhance the virulence of a pathogen or render a nonpathogen virulent 


Increase transmissibility of a pathogen 


Alter the host range of a pathogen 


Enable evasion of diagnostic/detection modalities 


Enable the weaponization of a biological agent or toxin 


Any other potentially harmful combination of experiments and agents 


Plants 


Seed stocks Report on the source of all seed stocks or other plant material used. If applicable, state the seed stock centre and catalogue number. If 
plant specimens were collected from the field, describe the collection location, date and sampling procedures. 


Novel plant genotypes Describe the methods by which all novel plant genotypes were produced. This includes those generated by transgenic approaches, 
gene editing, chemical/radiation-based mutagenesis and hybridization. For transgenic lines, describe the transformation method, the 
number of independent lines analyzed and the generation upon which experiments were performed. For gene-edited lines, describe 
the editor used, the endogenous sequence targeted for editing, the targeting guide RNA sequence (if applicable) and how the editor 
was applied. 


Authentication Describe any authentication procedures for each seed stock used or novel genotype generated. Describe any experiments used to 


assess the effect of a mutation and, where applicable, how potential secondary effects (e.g. second site T-DNA insertions, mosiacism, 
off-target gene editing) were examined. 


ChIP-seq 


Data deposition 


Confirm that both raw and final processed data have been deposited in a public database such as GEO. 


Confirm that you have deposited or provided access to graph files (e.g. BED files) for the called peaks. 


Data access links For "Initial submission" or "Revised version" documents, provide reviewer access links. For your "Final submission" document, 
May remain private before publication. | provide a link to the deposited data. 


Files in database submission Provide a list of all files available in the database submission. 


Genome browser session Provide a link to an anonymized genome browser session for "Initial submission" and "Revised version" documents only, to 


(e.g. UCSC) enable peer review. Write "no longer applicable" for "Final submission" documents. = 
= 

Methodology 
© 

Replicates Describe the experimental replicates, specifying number, type and replicate agreement. E 
© 

Sequencing depth Describe the sequencing depth for each experiment, providing the total number of reads, uniquely mapped reads, length of reads and G 
whether they were paired- or single-end. a 

D 

Antibodies Describe the antibodies used for the ChIP-seq experiments; as applicable, provide supplier name, catalog number, clone name, and S 
lot number. = 

3 

© 

Peak calling parameters Specify the command line program and parameters used for read mapping and peak calling, including the ChIP, control and index files wn 
used. = 

3 

Data quality Describe the methods used to ensure data quality in full detail, including how many peaks are at FDR 5% and above 5-fold enrichment. 9% 
ka 

Software Describe the software used to collect and analyze the ChIP-seq data. For custom code that has been deposited into a community 


repository, provide accession details. 


Flow Cytometry 


Plots 


Confirm that: 


The axis labels state the marker and fluorochrome used (e.g. CD4-FITC). 


The axis scales are clearly visible. Include numbers along axes only for bottom left plot of group (a ‘group’ is an analysis of identical markers). 


All plots are contour plots with outliers or pseudocolor plots. 


A numerical value for number of cells or percentage (with statistics) is provided. 


Methodology 

Sample preparation Describe the sample preparation, detailing the biological source of the cells and any tissue processing steps used. 

Instrument Identify the instrument used for data collection, specifying make and model number. 

Software Describe the software used to collect and analyze the flow cytometry data. For custom code that has been deposited into a 
community repository, provide accession details. 

Cell population abundance Describe the abundance of the relevant cell populations within post-sort fractions, providing details on the purity of the 
samples and how it was determined. 

Gating strategy Describe the gating strategy used for all relevant experiments, specifying the preliminary FSC/SSC gates of the starting cell 


population, indicating where boundaries between "positive" and "negative" staining cell populations are defined. 


Tick this box to confirm that a figure exemplifying the gating strategy is provided in the Supplementary Information. 


Magnetic resonance imaging 


Experimental design 


Design type Indicate task or resting state; event-related or block design. 


Design specifications Specify the number of blocks, trials or experimental units per session and/or subject, and specify the length of each trial 
or block (if trials are blocked) and interval between trials. 


Behavioral performance measures _| State number and/or type of variables recorded (e.g. correct button press, response time) and what statistics were used 
to establish that the subjects were performing the task as expected (e.g. mean, range, and/or standard deviation across 
subjects). 


Acquisition 
Imaging type(s) 


Field strength 


Sequence & imaging parameters 


Area of acquisition 


Diffusion MRI 


Preprocessing 


Preprocessing software 


Normalization 


Normalization template 


Used 


Noise and artifact removal 


Volume censoring 


Specify: functional, structural, diffusion, perfusion. 
Specify in Tesla 


Specify the pulse sequence type (gradient echo, spin echo, etc.), imaging type (EPI, spiral, etc.), field of view, matrix size, 
slice thickness, orientation and TE/TR/flip angle. 


State whether a whole brain scan was used OR define the area of acquisition, describing how the region was determined. 


Not used 


Provide detail on software version and revision number and on specific parameters (model/functions, brain extraction, 
segmentation, smoothing kernel size, etc.). 


If data were normalized/standardized, describe the approach(es): specify linear or non-linear and define image types used for 
transformation OR indicate that data were not normalized and explain rationale for lack of normalization. 


Describe the template used for normalization/transformation, specifying subject space or group standardized space (e.g. 
original Talairach, MNI305, ICBM152) OR indicate that the data were not normalized. 


Describe your procedure(s) for artifact and structured noise removal, specifying motion parameters, tissue signals and 
physiological signals (heart rate, respiration). 


Define your software and/or method and criteria for volume censoring, and state the extent of such censoring. 


Statistical modeling & inference 


Model type and settings 


Effect(s) tested 


Specify type (mass univariate, multivariate, RSA, predictive, etc.) and describe essential details of the model at the first and 
second levels (e.g. fixed, random or mixed effects; drift or auto-correlation). 


Define precise effect in terms of the task or stimulus conditions instead of psychological concepts and indicate whether 
ANOVA or factorial designs were used. 


Specify type of analysis: |_| Whole brain [_] ROl-based  [_] Both 


Statistic type for inference 


(See Eklund et al. 2016) 
Correction 


Models & analysis 


Graph analysis 


n/a | Involved in the study 


Specify voxel-wise or cluster-wise and report all relevant parameters for cluster-wise methods. 


Describe the type of correction and how it is obtained for multiple comparisons (e.g. FWE, FDR, permutation or Monte Carlo). 


Functional and/or effective connectivity 


Multivariate modeling or predictive analysis 


Functional and/or effective con nectivity Report the measures of dependence used and the model details (e.g. Pearson correlation, partial correlation, 


Graph analysis 


mutual information). 


Report the dependent variable and connectivity measure, specifying weighted graph or binarized graph, 
subject- or group-level, and the global and/or node summaries used (e.g. clustering coefficient, efficiency, 
etc.). 


Multivariate modeling and predictive analysis | Specify independent variables, features extraction and dimension reduction, model, training and evaluation 


metrics. 
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